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ABSTRACT 

The galaxy stellar-to-halo mass relation (SHMR) is nearly time-independent for z < A. 
We therefore construct a time-independent SHMR model for central galaxies, wherein 
the in-situ star formation rate (SFR) is determined by the halo mass accretion rate 
(MAR), which we call Stellar-Halo Accretion Rate Coevolution (SHARC). We show 
that the ~ 0.3 dex dispersion of the halo MAR matches the observed dispersion of 
the SFR on the star-formation main sequence (MS). In the context of “bathtub”-type 
models of galaxy formation, SHARC leads to mass-dependent constraints on the re¬ 
lation between SFR and MAR. Despite its simplicity and the simplified treatment of 
mass growth from mergers, the SHARC model is likely to be a good approximation 
for central galaxies with M, = 10® — that are on the MS, representing most 

of the star formation in the Universe. SHARC predictions agree with observed SFRs 
for galaxies on the MS at low redshifts, agree fairly well at z ^ 4, but exceed obser¬ 
vations at z 4. Assuming that the interstellar gas mass is constant for each galaxy 
(the “equilibrium condition” in bathtub models), the SHARC model allows calcula¬ 
tion of net mass loading factors for inflowing and outflowing gas. With assumptions 
about preventive feedback based on simulations, SHARC allows calculation of galaxy 
metallicity evolution. If galaxy SFRs indeed track halo MARs, especially at low red- 
shifts, that may help explain the success of models linking galaxy properties to halos 
(including age-matching) and the similarities between two-halo galaxy conformity and 
halo mass accretion conformity. 

Key words: cosmology: theory - galaxies: halos - galaxies: evolution - methods: 
N-body simulations - methods: luminosity function, mass function 


1 INTRODUCTION 

In the cold dark matter paradigm, galaxies form in dark 
matter halos. As cosmological simulations such as M illen- 
nium llSpringel et al.ll2005l: iBovlan-Kolchin et al.ll2009l) and 

m 


bright and faint galaxies th at are in excellent agreement 


Bolshoi i Klvpin, Trujillo-Gomez fc PrimackI 20111 )' resolved 
dark matter halos increasingly well, it has been a major goal 
to use such simulations to improve our understanding of the 
connection between h alos and the galaxies that they host. 
Abun dance matching llKravtsov et al.l2004l : IVale fc Ostrikerj 
|2004^ - which in its most basic form is just rank ordering 
galaxies by their stellar mass and assigning them to 
halos ranked by mass or peak circular velocity - leads to 
predictions of the galaxy autocorrelation functions for both 
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20061: iKlvoin. Truiillo-Gomez & PrimackI 

2011; 

Rodriguez-Puebla. Drorv & Avila-Reesd 

2012: 

Reddick et al.l 2013. and references therein). Abun- 


dance matching taking into account galaxy star formation 
rates also allows calculation of the typical relationship 
between the mass of dark matter halos and the stellar 
mass of the hosted galaxies (e.g. Moster, N aab fc White! 


l2013l : iBehroozi. Wechsler fc Conroy FoiSlJl . The resulting 
stellar-to-halo mass relation (SHMR) is remarkably similar 
at all redshifts between 0 and 4. This is consistent with the 
assumption that the average virial star formation efficiency 
(the star formation rate divided by the halo baryon accretion 
rate) is only a function of ha lo mass and not redshift from 
z = 4 to the present epoch llBehroozi, Wechsler fc ConrovI 
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l2013al '). This motivates us to develop a simple model in 
which the mass accretion rate (MAR) of dark matter 
halos determines the star formation rate (SFR) of their 
host galaxies. We call this the Stellar-Halo Accretion Rate 
coevolution (SHARC) assumption. Note that it is also 
possible to develop different galaxy-halo coevolution models 
that could also satisfy the time-independent SHMR by 
correlating SFRs to other halo assembly properties (e.g., 
halo formation time) but SHARC is a particularly simple 
one. 

The mass accretion rate of dark matter halos de¬ 
pends on the precise definition of the halos, which 
has been called into question in several recent papers. 
iDiemer. More fc KravtsTO ll2013ll argued that much of the 
mass evolution of dark matter halos is an artifact caused 
by the changing radius of the halo, a phenomenon that 
they call “pseudo-evolution.” In this paper we define the 
radius of the halo as the radius Rvir that encloses an aver¬ 
age density Avirpm, where pm is the mean matter density of 
the universe HmPc pc is critical density, and the redshift- 
dependent virial ove rdensity Avir(z) is given by the spher¬ 
ical collapse model llBrvan &: Normanlll998l i. Other popu¬ 
lar definitions are ii200m and i?200c, corresponding to en¬ 
closed densites of 200pm and 200pc respectively. For all these 
definitions, the rapid drop in background density as z de¬ 
creases is the main cause of the increase in halo virial radius 
and therefore a main cause of the increase in the enclosed 
mass, while the dark matter distributio n in the inter ior of 
the halo hardly changes at low redshift (IPrada et al. 


I2OO6I: 


IPiemand. Kuhlen fc Madau 200^; Cuesta et al.l 200g ). Re¬ 
cently More. Diemer fc KravtsoJ ll2015l l proposed that the 


best physically-based definition of halo radius is the “splash¬ 
back radius” Rsp « 2i?200m, where there is typically a sharp 
drop in the density; using this definition, there is actually 
more halo mass increase than for Rvir, R 200 m, or R 200 C. 

What is actually relevant to star formation of the 
central galaxy in the halo is the amount of gas that en- 
ters the halo and eve ntually reaches its central regions. 

I Wetzel &: Nagail (l2014l i have used adaptive refinement tree 
(ART) hydrodynamic galaxy simulations with a best res¬ 
olution of about 0.5 kpc to show that infalling gas decou¬ 
ples from dark matter starting at about 2 R 200 m and roughly 
tracks the growth of M 2 oom. Thus, they argue that pseudo¬ 
evolution is relevant to the accretion of dark matter, but 
not to that of gas. Further evidence that the gas falling into 
the central regions of halos roughly tracks the halo mass 
accretion rate is provided bv IPekel et al.l ll2013ll . who used 
a suite of ART hydrodynamic zoom-in galaxy simulations 
with Mvir/IO^^MQ = 0.1 to 2 at 2 : = 2 and best resolution 
of 35 pc, and found that the gas inflow rate is proportional 
to the halo mass accretion rate, with about half the gas 
penetrating to O.lRvir at redshifts 2 = 4 to 2, and with the 
fraction increasing from 2 = 2 to 1. Analysis of a subse¬ 
quent suite of ART hydrodynamic zoom-in galaxy simula¬ 
tions with better resolution and feedback showed that these 
simulated star-forming galaxies grow stellar mass at the 
same rate as the specific halo mass increase (IZolotov et al.l 
I2OI5I. Tacch ella et ah, in prep.). A sim i lar re sult has been 
reported in iGonzalez-Samaniego et ^ (l2014l l for galaxies 
formed in halos of ARir = 2 — 3 x lO^^M©. 

Star-forming galaxies are known to show a tight de¬ 
pendence of SFR on stellar mass, which is known as the 


“main sequence” of galaxy formation dSalim et abl 200?! : 
iNoeske et al.llioOTl : lElbaz et al.l[2007l : IPaddi et ah 2007 ). in 
analogy with the tight dependence on stellar mass of the 
properties of stars on the stellar main sequence. Park mat¬ 
ter halos also have a mass accretion rate that is roughly 
proportional to their mass, and we show in m\ of this pa¬ 
per that the dispersion of the halo mass accretion rate at a 
given halo mass is 0.3 to 0.4 dex, similar to the dispersion of 
the SFR on the main sequence. It was this equality that orig¬ 
inally motivates us to examine more closely the connection 
between mass accretion and star formation. 

Our analysis is limited to the connection between 
distinct halos and central galaxies only. One reason 
is that subhalos lose mass via tidal stripping, result- 
ing in negative values of accretion rates (see, e.g., 
Ivan den Bosch. Tormen fc Giocolil I 2 OO 5 I I. Satellite galaxies 
are also affected in other ways by their proximity to central 
galaxies. Therefore, studying the connection between sub¬ 
halo mass accretion and satellite SFRs is beyond the scope 
of this paper. In order to focus on main-sequence galaxies, 
we also discuss mainly dark matter halos of masses 10^^ to 
10^^ Mq. 

In ^ we will derive the SHMR for all SPSS 
central galaxies based on t he lYang et al.l (l2012f) and 
iRodrfguez-Puebla et al.l (l2015fl . Presently available data al¬ 
lows this to be done only for 2 ~ 0. In this paper, we make 
the simple assumption that this SHMR is valid at all red- 
shifts. The SHARC assumption will allow us to deduce the 
SFR for every halo in the Bolshoi-Planck simulation. When 
we compare these predictions with observations we find that 
they are in pretty good agreement from 2 = 0 to ~ 4, both 
for the SFR and its dispersion. 

Star formation is regulated by a complex interaction 
between gas inflows and outflows. Models that describe the 
basic processes that govern gas inflows and outflows and star 
formati on in galaxies are called “bathtub” models in the lit¬ 
eratu re (|Bgurdi4_et_^|201Cj;|Paye^J)inIntpr_&A3p]3enhei^^ 


2OIII. I2OI2I: iKriimholz, fc Pekell I2OI2I: IPekel et al.l 12013 


et al.l 20131 : iPekel fc MandelkCT 20141 : iForbes et alJ 
2 OI 4 I : Feldmarmll2015l : Mitra. Pave fc Finlatorl 2015ll . Such 
models consider both “equilibrium” conditions when the 
amount of gas in the interstellar medium (the “bathtub”) 
is constant because star formation equals net inflow, and 
non-equilibrium situations when the bathtub is filling or 
emptying. In the even simpler model in this paper we as¬ 
sume equilibrium at all times. We call this the Equilibrium 
and SHARC model, or E+SHARC. As is shown in ^ the 
time-independent SHMR is compatible with the equilibrium 
condition. 

This paper is organized as follows: ij2]describes the dark 
matter simulation used here and how we connect central 
galaxies to host halos. There we also make the simplifying 
assumption that the stellar-to-halo mass relation for central 
galaxies on the main sequences is independent of redshift 
and show how this allows us to infer SFRs from halo mass 
accretion rates, i.e., the SHARC assumption. In ^we ex¬ 
plore a very simple bathtub model, which we assume for sim¬ 
plicity to be in equilibrium at all times - i.e., the gas mass 
is constant - in order to identify and understand the condi¬ 
tions that are satisfied by the equilibrium time-independent 
SHMR model, i.e., the E-I-SHARC assumption. We show 
how net gas infall is connected to preventive feedback. As- 
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Figure 1. Halo mass accretion rates (MARs) from 2 = 0 to 2 = 3, from the Bolshoi-Planck simulation. The instantaneous rate is shown 
in black, and the dynamically time averaged rate in red. The gray band is the Icr (68%) range of the instantaneous mass accretion rates. 
All the slopes are approximately the same ~ 1.1 both for Mvir and Mvir.dyn- 


suming a power-law preventive feedback for Mvir ~ 10 ^^ Mq 
( representing virial shock heating of in-falling gas and the 
effects of AGN), we deduce mass-loading factors and their 
dispersion as a function of halo mass and redshift. In S|4]we 
deduce the SFRs and their dispersion implied by our model, 
and compare with observations of the SFRs and dispersion 
on the star-forming main sequence. Not surprisingly, we find 
that our simple model does not correctly predict the SFR at 
high redshifts 2^4, showing that the SHMR should change 
above 2 ~ 4. However, from 2 = 4 to 0, the SFR predic¬ 
tions from the time-independent SHMR model are in better 
agreement with the observed SFRs on the main sequence, 
especially if we use the most recent observations. In Sj5] we 
compare the cosmic star formation rate density and galaxy 
stellar mass function with observations up to 2 ~ 6, again 
finding that the SHMR assumption is disfavored at high 
redshift. In ij6]we calculate the metallicity of the interstel¬ 
lar medium given by our E+SHARC model and compare 
with observations, yet again finding that this model fails at 
high redshift. Finally, m summarizes our conclusions, dis¬ 
cusses their implications, and describes ways to increase the 
generality of the simplified model considered here. 

We adopt cosmological parameters Ha = 0.693, Hm = 
0.307, Hbar = 0.048, h = 0.678, Ug = 0.96, and as = 
0.829, consi stent with recent results from the Plan ck Col¬ 
laboration jPlanck Collaboration et al.ll20l3 . l2015l ). These 
are the parameters used in the Bolshoi-Planck simulation 
llKlvpin et al.|[^014l . Rodriguez-Puebla et al. 2015 in prep.), 
on which our results here are based; as noted above, our 
halo masses are defined using the spherical overdensity crite- 


Table 1. List of acronyms used in this paper. 


ART 

Adaptive Refinement Tree (simulation code) 

CSFR 

Cosmic Star Formation Rate 

IMF 

Initial Mass Function 

ISM 

Interstellar Medium 

GSMF 

Galaxy Stellar Mass Function 

MAR 

Mass Accretion Rate, Mvir 

SHARC 

Stellar-Halo Accretion Rate Coevolution 

E-I-SHARC 

Equilibrium+SHARC 

SDSS 

Sloan Digital Sky Survey 

SFR 

Star Formation Rate 

SHMR 

Stellar-to-Halo Mass Relation 

sMAR 

Specific Mass Accretion Rate, Mvir/Mvir 

sSFR 

specific Star Formation Rate, SFR/M* 


rion of iBrvan fc Norm"^ lll998l) . We also assume a lChabri^ 
ll2003l) IMF. Finally, Table [T] lists all the acronyms used in 
this paper. 


2 STELLAR-HALO ACCRETION RATE 
COEVOLUTION (SHARC) 

2.1 The Simulation 

We generate our mock galaxy c atalogs based on t he N- 
body Bolshoi-Planck simulation jKlvpin et al.l 12014) . The 
Bolshoi-Planck simulation is based on the ACDM cosmol¬ 
ogy with parameters cons istent with the latest results from 
the Planck Collaboration (IPlanck Collaboration et al.ll2015h 
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and run using the Adaptive Refinement Tree code (ART 
Kjavtsov, Klvpin fc KhokhlovI 1 19971 : iGottloeber fc Klvninl i 
2008h . The Bolshoi-Planck simulation has a volume of 
(250h“^Mpc)® and contains 2048® particles of mass 1.9 x 10® ^ 
Mq. Halos/subhalos and their merger trees were calculated 
with the phase-space temporal halo finder ROCK STAR « 
dBehroozi et al.l2013l : lBehroozi. Wechsler fc Wul2013l '). Halo 5 
masses were defined using spherical overdensities accord¬ 
ing to the redshift-dependent virial o verdensity Ayjrf^:) -1 

given by the spherical collapse model jBrvan fc Normaa 
Il998l ~). with Avir = 178 for large 2 and Avir = 333 ^ 
at 2 = 0 with our Hm- Like the Bolshoi simulation ^ ^ 

dKIypin, Truiillo-Gomez fc PrimackI I 2 OIII L Bolshoi-Planck | 
is complete down to halos of maximum circular velocity '^- 0.5 
^max - 55 km/s. 

In this paper, we calculate instantaneous halo mass ac- 
cretion rates from the Bolshoi-Planck simulation, as well as 
halo mass accretion rates averaged over the dynamical time 



log M„r [Mol log M.Lr [“ol log M.,,. [Mg] 


(A/vir,dyn), defined as 

/ dMvir \ _ Mvii.(t) - ALvir(f - fdyn) 

\ dt I dyn tuyn 


Figure 2. Scatter of halo mass accretion rates from 2 = 0 to 
2 = 3 from the Bolshoi-Planck simulation. As in Fig. [T] scatter 
(1) for the instantaneous rate is shown in black, and that for the 

dynamically time averaged rate in red. 


The dynamical time of the halo is tdyn(-2) = 
[GAvir(2:)pm]~^^^, which is 20% of the Hubble time. 
Simulations fe.g.. [Dekel et alll2009l '> suggest that most star 
formation results from cold gas flowing inward at about the 
virial velocity - i.e., roughly a dynamical time after the gas 
enters. As instantaneous accretion r ates for distinc t halo s 
near clusters can also be negative dBehroozi et al.l l2014h , 
using time-averaged accretion rates allows galaxies in these 
halos to continue forming stars. 

Figure [T] shows the instantaneous and the dynamical- 
time-averaged halo mass accretion rates as a function of 
halo mass and redshift, and Figure [2] shows their respective 
scatters. Even before converting halo accretion rates into 
star formation rates [ 112.311 . it is evident that both the slope 
and dispersion in halo mass accretion rates are already very 
similar to that of galaxy star formation rates on the main 
sequence. 


2.2 Connecting Galaxies to Halos 


The abundance matching technique is a simple and powerful 
statistical approach to connecting galaxies to halos. In its 
most sinmle form, the cumulative halo and subhalo mass 
functiorU and the cumulative galaxy stellar mass function 
(GSMF) are matched in order to determine the mass 
relation between halos and galaxies. In order to assign 
galaxies to halos in the Bolshoi-Planck simulation, in this 
paper we use a more general procedure for abundance 
matching. Recent studies have shown that the mean stellar- 
to-halo mass relations (SHMR) of central and satellite 
galaxies are slightly different, especially at lower masses 
where satellites tend to have more stellar mass compared 
to centrals of the same halo mass (for a more general 


discu ssio n see Rodrfguez;Pueblaj_^ror^^_^vil^Ree3 


201^;__ Rodrfeuez;Pueblaj_^vila^Reese^ DrorvI l2013l: 


Reddick et al.l l2013l : IWatson fc GonrovI l2013l : IWetzel et alJ 
2013h . Since we are interested in studying the connection 


^ Typically defined at the time of subhalo accretion. 


between halo mass accretion and star formation in central 
galaxies, for our analysis we derive the SHMR for central 
galaxies only. 

We model the GSMF of central galaxies by defining 
P(M*|Mvir) as the probability distribution function that a 
distinct halo of mass Mvir hosts a central galaxy of stellar 
mass M*. Then the GSMF for central galaxies as a function 
of stellar mass is given by 

POO 

<[>*,cen(M*) = / P(M.|Mvir)<[>h(Mvir)dMvir. ( 2 ) 

Jo 


Here, <[>h(Afvir) is the halo mass function and P(M*|Mvir) 
is a log-normal distribution assumed to have a scatter of 
ac = 0.15 dex independent of halo mass. Such a value 
is supported by the analysis of large group catalog s 


dVang. Mo fc van den Bosch l2009l : iReddick et ah _ 2013|L 


studi es of the kinematics of satellite galaxies ( More et al.l 
ImS), as well as cluster i ng an al ysis of large samples of 
galax ies dShankar et al.l I20l4 iRodrfguez-Puebla et al.l 
I 2 OI 5 II . Note that this scatter, ac, consists of an intrin¬ 
sic component and a measurement error component. 
At 2 = 0, most of the scatter appears to be intrin¬ 


sic, but 
redshifts 


that beco mes less and less true at higher 
(see, e.g.. iBehroozi. Conroy fc Wechsleil l2010l: 


2013bl: Leauthaud et ahl 

we do not deconvolve to 
as most of the observations 
these errors in their 


BehTqozij_^echsler&_Conrpy 

2 OI 2 I: [Tinker et al.ll20lfr Here 


remove measurement error 
that we will compare to include 
measurements. 

As regards the G SMF of central galaxies, we he re use 
the results reported in iRodrfguez-Pu ebla et al.l d2015 lL In a 
recent analysis of the SDSS DR7, Rodrfguez-Puebla et al.1 
d2015ll derived the total, central, and satellite GSMF for 
stellar masses from = IO^Mq to M* = IO'^^Mq based 
on the NYU-VAGC dBlanton et al.l l2005h and using the 
l/Knax estimator. The membership (central/satellite) for 
each galaxy was o btained from an updated ve rsion of the 
lYang et al.l d2007l l group catalog presented in lYang et al.1 
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Figure 3. Upper Panel: Stellar-to-halo mass relation (SHMR) 
for S PSS galaxies. The red curve is for all SPSS galaxies, 
from iBehroozi, Wechsler ConrovI ll2013lj^ abundance match¬ 
ing using the Bolshoi simulation. The black curve is for SPSS 
central galaxies, using the abundance matc hing method of 
iRodrfguez-Puebla. Avila-Reese &; Pro^ 112013^ applied to the 
Bolshoi-Planck simulation. The latter is what we use in the 
present paper, where we restrict attention to central galaxies. 
Bottom Panel: Halo-to-stellar mass relations. The dotted ver¬ 
tical line and the blue arrow indicate that galaxies below M* = 
IqIO.SMq are considered as main sequence galaxies, while some 
higher-mass galaxies are not on the main sequence. 


(|2012h . The corresponding SHMR is shown as the black 
curve in Figure [3l and the S HMR for all galaxies from 
iBehroozi, WechslCT fc ConrovI ((20133) is shown as the red 
curve. The difference between the two curves for halo 
masses lower than Mvir ~ IO^^Mq reflects the fact that 
the SHMR of centrals and satellite galaxies are slightly 
different as mentioned above. At halo masses higher than 
Mvir ~ IO^^Mq , this difference is primarily due to the dif- 
ferences between the GSM Fs used to derive these SHMRs, 
IBehroozi et al.l l2013l used iMoustakas et al] (l2013l i . When 
comp aring both GSMFs we f ind th at the high mass-end 
from iRodrfguez-Pueb la et al. (|2015|) is signih cantly differ¬ 
ent to the one deriv£hijMoustakas_et_alJ ((201^. In contrast, 
when comparing iRodrfguez-Puebla et ^ ( 2015l f GSMF with 


iBernardi et al.l 1 I 2 OIOII we h nd an excellent agreement, for a 
more general discussion see iRodrfguez-Puebla et al.l (l2015h . 
In less degree, we also find that the different values employed 
for the scatter of the SHMR explain these differences. 


2.3 Inferring Star Formation Rates From Halo 
Mass Accretion Rates 

A number of recent studies exploring the SHMR at dif¬ 
ferent redshifts hav e found that it evolves only slowly 
with time (see, e.g., LgmrLhnnd__et_al, 201^; iHudson et al.l 
I 2 OI 3 I : IBehroozi Wechsler fc ConrovI l2013bl . and references 
therein). For example, based on the observed evolution 
of the GSMF, the star fo rmation rate SFR, and the 
cosmic star formation rate, IBehroozi. Wechsler fc ConrovI 
(l2013bf l showed that this is the case at least up to 
z = 4 (cf. possibl e increased evolu t ion a t 2 > 4; 

Behrooz^^_Sil^ 201^ Fir]telsteh]_et al] l2015lf . Moreover, 


Behroozi. Wechsler fc ConrovI ( 2013ah showed that assum¬ 


ing a time-independent ratio of galaxy specific star forma¬ 
tion rate (sSFR) to host halo specific mass accretion rate 
(sMAR), defined as the star formation efficiency e, simply 
explains the cosmic star formation rate since 2 = 4. If we 
assume a time-independent SHMR, the star formation effi¬ 
ciency is the slope of the SHMR, 


M*/M* 


aiogM* 


Mvir/Mvir l91ogMv: 


(3) 


This equation simply relates galaxy SFRs to their host 
halo MARs without requiring knowledge of the underlying 
physics. (This is the main difference between the equilibrium 
solution we present below and previous “bathtub” models.) 
Our primary motivation here is to understand whether halo 
MARs are responsible for the mass and redshift dependence 
of the SFR main sequence and its scatter. Similar models 
have been explored in the past for different purposes, includ - 
ing generating mock catalogs (iTaghizadeh-Popp et al.ll2015fl 
and understanding th e different clu stering of quenched and 
star-forming galaxies (lBeckerll2015h . 

Using halo MARs, we operationally infer galaxy SFRs 
as follows. Let M* = M*(Mvir(t), t) be the stellar mass of a 
central galaxy formed in a halo of mass Mvir(t) at time t. 
In a time-independent SHMR, the above reduces to M* = 
M*(Mvir(t)). From this relation the change of stellar mass 
in time is simply 


dM« 

dt 


= R 


91ogM« dMvi 
d log Mvir dt 


(4) 


where /* = M*/Mvir is the stellar-to-halo mass ratio. 
Equation implies stellar-halo accretion rate coevolution, 
SHARC. The left panel of Figure U shows the resulting 
stellar-to-halo mass ratio, /*, derived for SDSS central galax¬ 
ies (see Section [23J- Consistent with previous studies, we 
Hnd that /* has a maximum of 0.03 at Mvu ~ IO^^Mq, 
and it decreases at both higher and lower halo masses. The 
product /* X e = dM*/dMvir will be shown as the black 
curves in Figure [5] below. 

In the more general case M* = M*(Mvir(t), 2 ), equation 
@ generalizes to 


dM* _ aM«(Mvir(f),2) dMv 


avir dMt{Mvii{t), Z) dz 
dt dMv\T dt dz dt’ 


(5) 


© 0000 RAS, MNRAS 000, 000-000 
























































































6 




[Mo] [Mo] 


Figure 4. Left Panel: Stellar-to-ha lo mass ratio for SPSS galaxies. Righ t Panel: Star formation efficiency. As in Figure (3] the red 
curves are for all SDSS galaxies, from lBehroozi, Wechsler &: ConrovI ll2013bl l. and the black curves are for SDSS central galaxies only for 
Bolshoi-Planck simulation. 


where the first term is the contribution to the SFR from 
halo MAR and the second term is the change in the SHMR 
with redshift. Although in this paper we assume a constant 
SHMR, the formalism that we describe below applies to this 
more general case. 

The relation between stellar mass growth and observed 
star formation rate is given by 

SFR = M*/(l-R), (6) 


where R is the fraction of mass that is returned as gaseous 
material into the interstellar medium, ISM, from stellar 
winds and short lived stars. In other words, 1 — R is the frac¬ 
tion of the change in stellar mass that is kept as long-lived. 
Here we make the instantan eous recycling approximation, 
with R = R(t) as derived in iBehro ozi. Wechsler fc ConrovI 
ll2013bl . §2.3) and consistent with the lChabrieiT i 2003l ] IMF. 
(In our model, for simplicity we take t to be the cosmic time 
since the Big Bang.) 


2.4 Star Formation Efficiency 


Note that the star formation efficiency e, equation ©, also 
quantifies galaxy stellar versus halo mass growth. The right 
panel in Figure |4] shows the star formation efficiency as a 
function of halo mass. Several features in this figure are 
worth discussing. As has been long established, the star 
formation efficiency decreases significantly from ARir ~ 
10 ^^Mq to Mvir ~ IO^'^Mq, implying strong differences be¬ 
tween galaxy and halo mass growth. Low-mass halos gain 
mass more slowly than low-mass galaxies. For high-mass 
halos, this trend is inverted: high-mass halos grow faster 
than high- mass galaxies. This is commonly called “down¬ 
sizing” (see Fonteno^et^ 20091 : IConrov fc Wechsl^l2009l : 
iFirmani fc Avila-Reesd l20inl . and references therein). Sec¬ 
ondly, Milky-Way sized halos, ARir ~ IO^^Mq, have a star 
formation efficiency of e ~ 1. Note that, for e = 1, galaxy 


mass growth becomes linearly proportional to the host halo’s 
mass growth (A/* oc A/vir). 

It is useful to rewrite the SFR as a function of the star 
formation efficiency, 

SFR = /. X e X ALvir/(l - R). (7) 


This robust new result of the SHARC assumption can also 
be written as 


R(ALviO 


SFR 
R ARir 


/* X e 

Hi-RY 


( 8 ) 


which can be termed the “virial star-formation efficiency.” 
Here the universal baryon fraction is defined as fh = 
flbar/flM, aud has a value of R = 0.156 with the Planck 
cosmological parameters adopted for this paper. 

While in the analysis described above SFRs are based 
on instantaneous halo mass accretion rates A/vir, we also de¬ 
rive relations when using mass accretion rates averaged over 
a dynamical time scale, ARir.dyn- Specifically, in equation 
© we substitute ARir by ARir^dyn, given by equation (1). 

We only expect equations © and m to apply to star¬ 
forming galaxies on the main sequence. For our SDSS cal¬ 
ibration sample this includes galaxies with A/* = 10® — 
IOIO.SMq, go it is this mass range, shown by the dotted line 
a blue arrow in Figure[3l where we use the SHARC assump¬ 
tion in the rest of this paper. Note that above these masses 
quenched galaxies are detached from the mass accretion-star 
formation correlation. 


2.5 Impact of Mergers 

Both in-situ star formation and galaxy mergers can con¬ 
tribute to the stellar mass growth of galaxies. But most 
mergers of galaxies with ARir > IO'^^A/q at low redshift 
are dry mergers, so in inferring the stellar mass growth we 
should not include halo mass growth due to dry mergers. 
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In addition to stars, merging halos may also contain dif¬ 
fuse circumgalactic medium. This will only contribute to 
the growth of the total baryonic content of the halo but not 
to the stellar mass growth of the central galaxy. We do this 
by multiplying the total stellar mass growth, M*, inferred 
naively from the halo mass growth, by the fraction of stellar 
mass growth of central galaxies that comes from star forma¬ 
tion, /sFR. Specihcally, we calculate in-situ star formation 
in central galaxies only as 

SFR= ^^/sFR(Mvn). (9) 

We parameterize .fs FR(Afvir) as a function of ha l o mass 
and redshift following iBehroozi. Wechsler fc ConrovI (l2013bl . 
equations 19-21 and Figure 10). 


3 SHARC + BATHTUB: EQUILIBRIUM 
ASSUMPTION, INFLOWS, AND 
OUTFLOWS 

The most important result from the previous section is that 
galaxy SFRs can be derived directly from halo MARs if the 
SHMR is time-independent, i.e., the SHARC model. In this 
section, we explore a very simple gas-regulated model in 
order to identify and understand the conditions that are 
satisfied by the SHARC model. As we will show later, the 
time-independent SHMR is compatible with the equilibrium 
condition where galaxies are regulated only between inflow, 
outflow and star-formation. 


Thus, the equation that governs the gas mass growth is 

Afgas.ISM = Mgrav.ISM+Afr.ISM—(1 —R)SFR— Mout,ISM.(10) 

It is more convenient to rewrite Mgas.iSM as a function of 
Mvir and SFR. To do so, we define £bB as the efficiency 
with which the inflowing cosmological baryons will penetrate 
down to the galaxy’s ISM, ?7w,ism as the mass loading factor 
of gas outflows, and ?7r,iSM as the mass loading factor of gas 
re-infalling. Hence, 

Afgas.ISM = /bfeffMvir “ (»?w,ISM “ 77r,ISM + 1 — R)SFR. (11) 

3.2 Equilibrium Condition 

In the equilibrium solution, galaxies are regulated only be¬ 
tween inflow, outflow and star-formation. The net change of 
the gas mass within the ISM is zero, Mgas.iSM = 0. Within 
this assumption, the SFR is given by 

SFR = ^TT^^vir. (12) 

77 -I- 1 — it 

Here we define rj = ryw.iSM — 77r,iSM as the net mass loading 
factor. 

At this point, it is worth mentioning why the approach 
followed in this paper is particularly relevant. This similar¬ 
ity between equations 0 and ini) is not a coincidence. It 
reflects the fact that a time-independent SHMR is compat¬ 
ible with the equilibrium condition, as we will show in the 
next section. 


3.1 Gas Equation 


We begin by defining the gas equation that describes 
the change of cold gas mass in the interstellar medium 
(ISM) of a galaxy. The model described in the fol¬ 
lowing section is a simple version of previous models 
discusse d in the literature, sometimes called “bathtub” 
mode l s (|Bouche_A-alJ|20ld; DaveiFmIatOT_&_0££enhuhn^ 


2 OIII. I 2 OI 2 I: iKriimholz fc Dekell I 2 OI 2 I: iDekel et al.l I 2 OI 3 I: 


Lilly et al. 2011^ Dekel fc Mandelkeil 2014: Forbes et al.l 


2 OI 4 I : Feldmarml l2015l : Mitra. Dave fc FinlatoJ 2015l l d^ 

scribing the basic processes that govern the different com¬ 
ponents in galaxies. 

The change of the total cold gas mass within the ISM 
of a galaxy, Mgas.iSM, is the result of the following physical 
mechanisms: 


(i) The rate at which the cosmological baryonic inflow 
material will reach the ISM of the galaxy, Mgrav.iSM. 
This process is assumed to be related to the gravitational 
structure formation of the halo, i.e., proportional to its 
mass accretion rate, M^vir. 


(ii) The rate at which the gas that was previously ejected 
due to outflows is re-infalling into the galaxy’s ISM, Ah, ism- 

(hi) The gas mass that is lost due to star formation 
corrected by the fraction of the material that is instanta¬ 
neously returned into the ISM, (1 — R) x SFR. 

(iv) The gas mass loss of the galaxy’s ISM ejected due to 
outflows, Mout.ISM. 


3.3 Inflows, Outflows and Re-infall 


Now we combine the equilibrium condition with the time- 
independent SHARC assumption, and call the combination 
E-bSHARC. That is, we combine equations (O and (1121) to 
give 

fhSeS f*e 

n + i-Ri-R' ^ 

It is illuminating how the above equation explicitly 
relates the parameters from the equilibrium condi¬ 
tion (left-hand side) to the SHARC assumption (right- 
hand side). Substantial progress has been made in 
mod eling the lef t han d side of equation (1131). see 
e.g.. iBo uche et ahl ll20 10ll: Dave. Finla tor fc Oppcnheimeij 


(l2f)l]}l : lMitra. Dave fc Finlatorl ( 2015l 'l. Nonetheless it is still 

challenging mainly because it involves many physical pro¬ 
cesses that are poorly constrained. In the time-independent 
SHMR model, however, the value of the left-hand side in 
equation m is constrained by the known value of the 
right-hand side. While the above does not give any infor¬ 
mation of the halo mass and redshift dependence of £e« and 
rj separately, this is possible if one uses prior information 
based on models of galaxy formation. For example, out¬ 
flows and re-infall are thought to be more relevant to halos 
of mass Mvir ~ 10 ^^Mq than in halo more massive than 
Mvir ~ IO^^Mq (as we will see later in Figure |6}. If we as¬ 
sume that such galaxies accrete at the maximum rate, i.e., 
£eS ~ 1, a crude estimation for the net mass-loading factor 
is given by 77 oc (/*e)“^. For lower-mass galaxies, preventive 
feedback may lead to £b« << 1, as we mention below. The 
situation is again different at higher masses, where the ac¬ 
cretion of cold gas is diminished and outflows and re-infall 
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Figure 5. Effective penetration parameter at 2 = 0,1, 4, and 6. The vertical axis is either, /* x e, /t, x S^, fb x Sq or fb x £’eff - Here 
1 — is the fraction of gas that never makes it into the galaxy as a result of “preventive feedbacks” associated with gas heating by 
virial shocks equation I15II and AGN feedback (5q, equation 1 1611 . As explained in the text, is based on results of other papers. 


are less relevant, i.e., ~ 0. Thus at high mass we expect 

that fbSeS = f*e. 

The penetration parameter £es is the result of various 
forms of preventive feedback, including: 

(i) Photoionization heating, fph. This term only affects 
very low-mass halos, Mvir ~ IO^Mq, which we do not 
consider in this paper. 

(ii) Heating of the inflowing cosmological baryons due to 
energetic winds, Winds are more significant in halos 
lower than Mvir ~ 10 ^^ Mq. 

(iii) Heating of inflowing gas as it crosses virial shocks, 
£h_- This term becomes relevant in halos more massive than 
Mvir ~ 10^^ 

(iv) Any process in massive halos that prevents cool¬ 
ing flows from reaching the central galaxy’s ISM, for ex¬ 


ample because of maintenance-mode feedback from super- 
massive black holes {£q), which also keeps quenched galaxies 
quenched. This term becomes more relevant in halos more 
massive than Mvu ~ 10^^. 

Following IPave. Finlator fc Qppenheimeil (l2012l l the result¬ 
ing £eff is given by 

~ £ph ^ £vi X £h ^ (14) 

For simplicity, we will ignore the impact of £ph and - i.e., 
we assume that £pb ~ fw ~ 1- This is well justihed given 
the halo mass scales Mvir ^ 10 ^^Mq analyzed in this paper. 

We now describe the functional forms we assume for 
£b and gg. From analysis of hydr odynamic simulations, 
iFaucher-Giguere. Keres fc Mai (1201 il l derived the halo mass 
and redshift dependence of £b given by 

a(Mvir, = min 11, 0.47 (l^^) } 
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Figure 6. Left Panel: Net mass loading factor, 77 = r^wjSM ~ "^rjSM a function of halo mass at 2 = 0,1, 4, and 6 , obtained assuming 
preventive feedback described by = 5h x £’q. The cal culated dispersion is shown at 2 : = 0 and 6 . We compare our result s with 
observational constraints on t he mass loading f actors from iBouche et alj <|2012|^ for a sample at 2 : ~ 0.1, iKacprzak et al.l <|2014|') for a 
galaxy sample at 2 : 0.2 and ISchroetter et alJ l|2015l'l for galaxies at 2 : 0 . 8 . Note that these observational constraints are referred to 

outflowing mass loading factors. Right Panel: Net mass loading factor and its dispersion as a function of galaxy stellar m ass. Empirical 
constraints on the mass loading factor based on an analytical model for galaxy metallicity in |Lu. Blanc & B enson] l|2015f) are s hown as 
the magenta and orange lines when using gas phase metallicity constraints from iMaiolino et alj |j2008^ and IZahid eT^ ll2013h z ~ 0 , 
respectively. Similarly above, these results are referred to outflowing mass loading factors. 


(see also IPave. Finlator fc Oppenheimerl I 201 II . I2ni2h . Fig¬ 
ure [ 5 ] shows the feedback parameter fh multiplied by the 
universal baryon fraction, /b, as the blue long dashed-line 
at z = 0,1,4 and z = 6. In the same hgure, we show the 
halo mass dependence of /* x e as the black solid line. Recall 
that we assume that /* x e is independent of redshift. 

If /* X e = /b X feff, then the mass loading factor rj = 
0. The difference between these two quantities is therefore 
related to rj. The point of maximum approach between these 
two curves is when ft x e reaches its peak value at Milky 
Way sized halos, Mvir ~ 10^^ Mq, where /b x fh is a factor 
of ~ 1.2 higher than /* x e at z = 0. At z = 6 the situation is 
qualitatively different and /b x fh is a factor of ~ 2.5 higher 
than /* X e in Mvir ~ 10^^ Mq halos. Then the mass loading 
factor should increase at high redshift. 

Halo mass quenching is more relevant for high mass 
halos. This imposes the constraint that any functional form 
proposed for Sq should reproduce the fall off at higher masses 
of the term /*£. Given the uncertain redshift dependence of 
£q, we will assume for simplicity that it is independent of 
redshift. The functional form £q that describes the fall-off 
of /*£ at z = 0 is given by 

f,(Mvn)=min|l,0.85(^^) “'j . (16) 

Note that at z = 0 for halos more massive than ~ 10^^ Mq, 
£es ~ £ X ft/fb- Such a fall-off is thus necessary in or¬ 
der to make SHMR-|-equilibrium assumptions work, in other 
words, equation inj. The green long dashed-dotted lines in 
Figure [ 5 ] show £q. At higher redshifts £ea > £ x /*//b imply¬ 
ing that the mass-loading factor becomes more important at 
high redshifts in high mass galaxies. 


Next, in equation ca) we use the functional forms de¬ 
scribed in equations (na and (na to deduce a relation for 
the net mass loading factor: 


V ■ 


fb £eff(Mvir,z) 
/*(Mvir) £(Mvir) 


- 1 


(1-R). 


(17) 


The left hand panel of Figure [ 6 ] shows the net mass loading 
factor, 77 = » 7 w,iSM — TJr.iSM, as a function of halo mass at z = 
0,1, 4 and 6 . Note that the generic redshift evolution of 77 is 
governed by the evolution of feff- For halos less massive than 
~ IO^^ ^Mq, Figure [ 6 ] shows that the mass loading factor 
approximately scales as a power law with a power that is 
roughly independent of redshift, 77 oc . Equivalently, 

we And that for galaxies with stellar mass below ~ 10 ®'^Mq 
the mass loading factor scales as 77 oc Mass loading 

factors are predicted to be very small for halos more massive 
than 10^^ Mq, especially at low redshifts. For comparison 
we include observational constra ints on the mass loading 
factors from iBo uche et ^ ll2012|] for a sample at z ~ 0.1, 
Kac 2 rzak_et^ I ll2014ri for a galaxy sample at z ~ 0.2 and 
Schroetter et al.l 1 2015l] for galaxies at z ~ 0.8. Empirical 


constraints on the mass loading factor based on an analytica l 
model for galaxy metallicity in ILu. Blanc fc BensorJ (l2015tl 
are shown as the magenta and o range lines when using gas 
phase metal li city c onstraints from iMaiolino et al.l ll2008l] and 
IZahid et al.l (l2013l] at z ~ 0, respectively. Note that these 
comparisons are referred to outflowing mass loading factors. 
Nevertheless, at lower masses, this comparison is fair since 
most of the outflowing gas is the most relevant contribution 
to the net mass loading factor. 

In this Section we presented a simple framework that 
clarihes how the net mass loading factor is connected to 
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preventive feedback in the context of the equilibrium time- 
independent SHMR model. As long as the SFR is driven 
by MAR these assumptions can be generalized in the same 
framework, as we mention briefly in the discussion section. 


4 SPECIFIC STAR FORMATION RATES 
FROM SHARC 


4.1 SHARC Compared with Observations 


We have now collected together all the tools needed to fol¬ 
low several aspects of galaxy evolution while galaxy stel¬ 
lar masses are in the range M* = 10®Mq to 
We start by showing the evolution in the slope and zero- 
point of the star-forming main sequence inferred by the 
time-independent SHMR (SHARC model) in Figure [T] Re¬ 
call that when assuming a time-independent SHMR, stellar 
mass growth can be inferred directly from halo mass accre¬ 
tion rates via M* = /* x e x Mvir, with the corresponding 
SFR = M*/(l — R). Black solid lines show results using 
instantaneous mass accretion rates, Mvir, in equation 
Red solid lines show the SFRs when using mass accretion 
rates smoothed over a dynamical time scale, Mvir.dyn, in¬ 
stead. The gray band indicates the intrinsic scatter around 
the star-forming main sequence when using Mvir. Note that 
our model sSFRs were corrected in order to take into ac¬ 
count the contribution of mergers to stellar mass growth, 
as explained in 323] We show the resulting sSFRs without 
this merger correction with the black and red dashed lines 
when using Mvir and Mvir.dyn respectively. Note that the 
contribution from mergers becomes more important for red- 
shifts 2 < 0.5. Hereafter, we will focus our discussion on 
the merger-corrected results, also shown as the solid lines in 
Figure [Hi 

Both Mvir and Mvir.dyn produce similar relations at all 
redshifts. This is expected due to the similarities shown be¬ 
tween Mvir and Mvir.dyn iu Figure [T] We note that the re¬ 
sulting slope of model the star-forming main sequence when 
using Mvir at 2 = 0 is 0.73 and increases as a function of red- 
shift to a value of 0.87 at 2 = 6. (Similar slopes are derived 
for Mvir.dyn.) TMs is Consiste nt with observed s l opes derived 
from SPSS galaxies (s ee e.g, lEIbaz et al.ll200^ : IZahid et al.l 
l2012l:ISalim' .l200'<tl as well as from high-redshift galaxies 
(see e .g, Santini et al. I I 2 OO 9 I : iKarim et al.lf20lil : iReddv et al.l 

l2012lb 


In Figure [Tj w e reproduce the best fit reported in 
ISneagle et al.l (|2014l ') to the star formation main sequence 
as the orange curve , as well as more recent observations. 
ISpeagle et al.l ll20l4 ) used a compilation of 25 observations 
from the literature to study the star formation main se¬ 
quence from 2 = 0 to 6. The authors carefully calibrated 
observational SFRs, correcting for different assumptions re¬ 
garding the stellar IMF, SFR indicators, SPS models, dust 
extinction, emission lines and cosmology, among the most 
important calibrations. Hence, their best-fitting model rep¬ 
resents a robust inference of the redshift evolution of the 
star-forming main sequence. Our deriv ed star-forming main 
sequences are in good agreement with lSneagle et ^ (l2014l i 
and within the Icr intrinsic scatter of our relations at almost 
all redshifts. Note, however, that there are some system¬ 
atic deviations between our model predictions and the ob- 
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Figure 9. Scatter of the sSFR for main-sequence galaxies pre¬ 
dicted in our model. 


servations. While these differences could be due to redshift- 
dependent systematic biases in the observationally-inferred 
sSFRs, it is interesting to discuss these differences in the 
light of the constant SHMR model. 

First, the observed sSFRs of galaxies at 2 > 4 are sys¬ 
tematically lower than the time independent SHMR model 
predictions. These differences increase at 2 = 6. The dis¬ 
agreement between the constant SHMR predicted SFRs and 
the observations implies that the changing SHMR must be 
used, as in equation m, at least at high redshift. 

Between 2 = 4 and 2 = 3 the observed star-forming 
sequence is consistent with the SHARC predictions. Between 
2 = 2 and 2 = 0.5, the observed sSFRs are slightly above 
the SHARC predictions. This departure occurs at the time 

of the peak value of the cosmic star formatio n rate. _ 

_ After the compilation carried out by ISneagle et al.l 

||2014D . new determinations of the sSFR have been pub¬ 
lished, particularly for redshifts 2 < 2.5. In FiguresEi 
we reproduc e new d ata published in |WhitaImr_et_^y_| 


,_, ,_, _ 2OII; 

lllbert et~^ (l2015fl: ISalmon et al] ( 2015fl : Schreiber et al.l 
( 2015h arid Tasca et al. ( 2015h . This new set of data agrees 
better with our model between 2 = 2 and 2 = 0.5, im¬ 
plying that the time-independent SHMR (SHARC assump¬ 
tion) may be nearly valid across the wide redshift range from 
2 ~ 4 to 2 0, a remarkable result. However, it is not clear 

whether this is valid si nce the newer observ ations have not 
been recalibrated as in lSneagle et al.l ll2014h . 


4.2 Scatter of the sSFR Main Sequence 

We now turn our discussion to the scatter of the star-forming 
main sequence, displayed in Figure [21 When using Mvir, the 
scatter is nearly independent of redshift and it increases 
very slowly with mass for 2 < 2. The value of the scat¬ 
ter is not symmetric and has values of ct ~ 0.35 — 0.45 dex. 
In contrast, for 2^2 the scatter decreases with mass. In¬ 
stead using Mvir.dyn produces a scatter more symmetric and 
practically independent of mass and redshift below 2 = 2. 
The scatter takes a value of cr ~ 0.35 dex. At high red¬ 
shifts, the scatter decreases with increasing mass. We cannot 
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Figure 7. Redshift evolution of sSFRs derived in the time-independent SHMR model (SHARC assumption). The red and black curves 
are the medians of the sSFRs from the dynamically-time-averaged and instantaneous mass accretion rates, respectively, the dispersion 
about the median sSFRs with the gray band calculated from instantaneous mass accretion rates . . The se are corrected for mergers (see 
while the corresponding dashed curves are not. These are compared with the lSpeagle et al.l ||2014|~I summary o f observed sSFRs on 
the main sequence (orange curve), and also with other recent measurements. Note that the lwMta^r et all ll2ni4l i measurements from 
z = 0.5 — 1 . 0, list e d as Whitaker+14, a re shown in bot h z = 0.5 and z = 1 pan els. I n addition, we incl ude star-forming data from 
lllbert et all ll2015l l ; ISalmon et all ll201.'jl ; ISchreiber et al.l ll2015l l ; iTasca et all ll2015l 'l and lR,eddv et al.l 120121 ). 


make a direct comparison with observations due the uncer¬ 
tainties affecting the measurements of both stellar masses 
and SFRs. Nevertheless, attempts to deconvolve the intrin¬ 
sic scatter from measurement errors, particularly for SDSS 
galaxies, have found that the star-formi ng main sequence 
has a scatter of a 0.3 dex (see, e.g., Salim et al.l l2007l : 
ISpeagle et al.|[2014l : ISchreiber et al.l 12015 ). New results on 
star-formation rates from CANDELS optical-IR colors also 
support a main-sequence scatter of ~0.3 dex that is remark¬ 
ably uniform from z = 0.5 to 2.5 at all masses 10®Mq and 
above (Fang et ah, in prep.). The behavior of the SFR dis¬ 
persion can perh aps be under stood as reflecting the Central 
Value Theorem llKelsorJl20l4l applied to the halo MAR. 

In this paper, we are using the MARs for all distinct 
dark matter halos to predict the SFR on the galactic main 


sequence and its scatter. What if we instead only used the 
halos that host star-forming central galaxies? In work in 
progress, we have found that doin g this at z ~ 0 using age 
matching jHearin Sz Watsonll2013l ) and similar methods re¬ 
sults in a somewhat smaller scatter in the predicted sSFR 
of about 0.3 dex, in even better agreement with observa¬ 
tions. It is not clear whether this will also be true for z > 0, 
however. 
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Figure 8 . Specific star formation rates as a function of redshift 2 for stellar masses M* = 10®, 10®'®, 10^® and 1O^® ®M 0 from time- 
independent SHMR model. The red and black curves are the sSFRs, from both dynamically-time-averaged and instantaneous mass 
accretion r ates, respectively, wi th the gray band representing the dispersion in the latter. Both are corrected for mergers . The orange 
curve is thelSoea gle et ahl ||2014| ) summary of observ ed sS FRs on the main seq uence. Observations from lWhitaker et al.l l|2Q14lh[ilbert et alJ 
ll2015l^ , ISalmon et al.l ll2015lb Schreiber et alj <l2015h and iTasca et al.1 ll2015l) are also included. 


5 SHARC EVOLUTION OF THE COSMIC tion of the GSMF is consistent with observational inferences 

STAR FORMATION RATE & STELLAR at different redshifts. 

MASS FUNCTION 


The SHARC model can be used to calculate the total cos¬ 
mic star formation rate (CSFR) as a function of time. This is 
shown in Figure [TOl which plots the results of using M^ii but 
similar results are obtained if Mvir.dyn is used instead. For 
comparison, we also reproduc e a com pilation presented in 
iBehroozi. Wechsler fc ConrovI (l2013bll . including data from 
UV, UV-fIR, IR, Ha and 1.4 CHz, as well as a recent fit to 
observ ations including UV and IR by iMadau fc DickinsonI 
ll2014h . The peak of our CSFR occurs a t 2 ~ 3 which is ear¬ 
lier than th e fit of Madau fc DickinsonI ll2014ll and the data 
compiled in IBehroozi. Wechsler fc ConrovI 1 2013bfl . peaking 


at 2 2. 


Figure [To] shows that the CSFR from the SHARC model 
is higher than most of the observations at 2 > 2, although as 
we showed in the previous section, the SFR predicted by th e 
model agrees between z=3 and 4 with lSneagle et al] (l2014h . 
However, at 2 > 4, where the model predicts high SFRs, 
this is a failure of the constant SHMR condition. 

Figure [TT] shows the redshift evolution of the CSMF cal¬ 
culated from the time-independent SHMR, SHARC assump¬ 
tion. In the same Figure we compare to some observational 
inferences as indicated in the caption. The predicted evolu- 


6 SHARC+EQUILIBRIUM BATHTUB: 
METALLICITIES 

In §4 and §5, we only used the time-independent SHMR 
(SHARC). Now we return to the bathtub model §3. Using a 
time-independent SHMR plus the Equilibrium assumption, 
E-I-SHARC, we can predict the metallicities in the ISM. 
Metallicity is defined as Zism = Mz,isM/Mgas,iSM, where 
Afz,isM is the mass of metals in the gas phase within the 
ISM and Mgas.iSM the total cold gas mass. The change of 
the metal mass within the ISM of a galaxy is given by 

Mz.ISM = y X SFR -I- ZlGMfbSe(lMvir + 
(.^r.ISM^r.ISM — (1 — R)^ISM — ^w,ISM??w,ISM) SFR. (18) 

The yield y in equation m is the metal mass in the 
gas phase formed and returned to the ISM per unit SFR. 
We assume that the metal yield is instantaneous, and use 
2 //(l — R) = 0.054 as derived in iKrumholz fc Dekell (l2012fl 
for solar metallicites. The terms Zigm, Zj-jsm, and Z„jsm 
are the metallicities of the intergalactic medium, the re-infall 
of previously ejected material, and the inflows of the ISM, 
respectively. 

Let Or and a™ be defined as ^r,iSM = Or^iSM and 


© 0000 RAS, MNRAS 000, 000-000 




















































Stellar-Halo Accretion Rate Coevolution 


13 


CO 

I 

O 

a 


I 

u 


® 0 



1 +z 


Figure 10. Cosmic star-formation rate (CSFR) as a function of 2 , for both our model (with results integrati ng down to log M*/Mq = 
8 , 9, and 1 0 shown as magenta, black, and green curves), and for fits bv iBehroozi. Wechsler & : Conro'^ l|2013 b|i (shown as a long-da shed 
red curve). [Madau &; DickinsonI ll2014l) (dotted blue curve), and sets of compiled observations l^ehro^i. Wechsler Conrovl[2013bl) . As 
discussed in the text the predicted CSFR in the SHARCK are probably related to the failure of the time-independent assumption. 


•^w,iSM = 0(-wZisM‘ Also, let the metallicity of inflowing ma¬ 
terial from the intergalactic medium be some fraction qigm 
of the galaxy’s ISM (Ziqm = Q^igm^ism). If we use the fact 
that ZisM = Mz,isM/Mgas,iSM in equation (HHj, it follows 
that the metallicity in the ISM is 


•^ISM = 


y X SFR — ZisM X Mgas.isM 


(OmTiw.ISM — Q!T-»7r,ISM + 1 — -R)SFR — aiGM/b^’effMv 


(19) 


In what follows we will assume that the outflowing 
metallicity is equal to the ISM metallicity (a-uj ~ 1). Also, we 
will assume that the galaxy ISM metallicity changes slowly 
compared to the re-infall time, so that Or ~ 1, and that 
the metallicity in the IGM is close to zero, i.e., ctigm ~ 0. 
If the galaxy’s metallicity evolves only slowly with time 
(^iSM ~ 0 ), the above equation can then be re-written as 


•^ISM = 


y 

R' 


( 20 ) 


This familiar equation is simila r to that in 
iDave, Finlator fc Qppenheimeil ll 2012 li . Under these 
assumptions the metallicity of the ISM is controlled by 


the net mass-loading factor 77 , which is itself controlled by 
the efficiency £ea (see equation fT3l) . In other words, the 
evolution of the metallicity in a galaxy’s ISM is driven only 
by the efficiency at which the baryons penetrate down to 
the galaxy. Note that this is no longer valid if e n riched 
outflows are considered, i.e., a™ > 1 . iFeldmarmI (l2015fl 
studied the general case and concluded that even if galaxies 
are not in strict equilibrium the outflowing metallicity is 
close to that of the galaxy’s ISM. Using equation , we 
can thus write Zism as 


.^ISM = 


/* X £ y 

fh X £eS 1 — R 


( 21 ) 


Note that a direct consequence of the equilibrium con¬ 
dition is that the scatter in Zism is simply a consequence of 
the scatter in the SHMR, i.e., that halos of the same mass 
have a range of values of /*. This is an important conclu¬ 
sion because if the observed scatter in Zism is similar to the 
intrinsic scatter of the SHMR, that would provide further 
evidence for the equilibrium condition. 

The resulting metallicities are compared with observa¬ 
tions in Figure[T2] Agreement would support the E-I-SHARC 
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Figure 11. Redshift evolution of the galax y stellar mass function 
sumpt ion), and a s observationally-derived in iRodnguez-Puebla et al. 
12005) at 2 : = 1; iMarchesini et al] ||2009|) at 2 : = 2: iMortlod^^alT 
( 2013 ) at 2 : = 6. 


(GSMF) as predic ted by the time-indepe nden t SHMR (SHARC as- 
112^^|) at 2 : = 0. l :lMoustal^s et alj ll2013l) and lPerez-Gonzalez et ^ 
1120111 at 2 = 3; fLee~et'^n ~ l201^) at 2 = 4 and 5: and IStar^’eT^ 


assumption, our preventive feedback assumptions, and our 
simple metallicity treatment. The model predictions are ac¬ 
tually in good agreement with some of the observations from 
2 = 0 to ~ 2, but it is hard to draw strong conclusions be¬ 
cause of the disagreements between different observations. 
However it is clear that our constant SHMR equilibrium 
model predicts metallicities that are much higher than ob¬ 
served at 2 3.5, which is a consequence of the model’s 

overprediction of the SFR at high redshift. 

7 CONCLUSIONS 
7.1 Summary 

In the present paper we have investigated to what extent 
the mass accretion rate of the host halo controls the rate 
of star formation of galaxies on the main sequence of star 
formation. We were motivated by the realization that the 
halo mass accretion rate (MAR) and its scatter, shown 
in Figures [T] and [21 bear a remarkable resemblance to the 
star formation rate on the main sequence and its scatter. 
In order to connect these phenomena, we have considered 


an extremely simple - no doubt oversimplihed - approach 
in which we made the crucial assumption that the stellar- 
to-halo mass relation (SHMR) for central galaxies in dark 
matter halos, deduced from SDSS observations and shown 
in Figure (3] remains valid at all redshifts. We called this 
Stellar-Halo Accretion Rate Coevolution (SHARC) assump¬ 
tion. We showed in §2 that the SHARC assumption allows 
derivation of galaxy SFRs from halo MARs. This robust new 
result can also be expressed as equation for the “virial 
star formation efficiency”, which only depends on halo mass, 
Mvir- In [j3] we showed that the SHARC assumption, i.e., 
a time-independent SHMR, is compatible with the equilib¬ 
rium condition that determines the amount of gas reaching 
the interstellar medium due to preventive feedbacks to the 
net mass-loading factor. We call this E-I-SHARC. Assuming 
reasonable preventive feedbacks based on simulations allows 
calculation of the mass loading factors and their dispersions 
shown in Figure jS] 

The specific star formation rate (sSFR) and its disper¬ 
sion predicted by the SHARC assumption are shown in Fig¬ 
ures IZllHl and (9) Despite the simplicity of our assumptions, 
the resulting predictions are in rather good agreement with 
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3 12. Metallicit i es as a function of stellar ma s s, both as predicted by the E+SHARC model and as observed bvlTrem onti e t alj 
iMaiolino et alJ ll2Q08h:fMannucci et al.N2010l'l:IZahid et alJ ll2Q12h:lHenrv et al~ l|2013l 'l: rSanchez et alJ l|2013l'l : [^eidel et al.ll|20l3^ 
Troncos^^^a?Tll2Q14h : IWuvts e^^ ll2014l^ : IZahi^^t al.l ll2014l^ and lSand^rs^rmi2013) . Conversion fromZigM to 12+log(0/H) following 


Fi| 

ilioo. 


PdldmamTil^Ol^ . The E+SHARC model predictions are in good agreement with observations below z <2 but disagree at higher redshifts 


probably reflecting the failure of the equilibrium condition. 


observations from z = 0 out to z ~ 4, especially if the most 
recent observations are used. At a ^ 4 the predicted sSFRs 
are systematically higher than observations, implying that 
the redshift-independent SHMR assumption breaks down at 
z >4. 

Figures [TU] and ED show that the cosmic star forma¬ 
tion rate density (CSFR) and galaxy stellar mass function 
predicted by the SHARC assumption is in good agreement 
with observations at z ^ 1, but the CSFR density is higher 
than most observations at higher redshifts. This disagree¬ 
ment at higher redshifts again arises because the redshift- 
independent SHMR assumption is invalid at higher red¬ 
shifts. 

Now, assuming SHARC plus an equilibrium bathtub 
model (E+SHARC) and a few additional plausible assump¬ 
tions regarding preventive feedback and inflow and outflow 
metallicites, the model also predicts ISM metallicities as a 
function of galaxy stellar mass and redshift. These predic¬ 
tions are compared with observations in Figure ED The pre¬ 
dictions are in good agreement with at least some of ob¬ 
servations z ^ 2, although the scatter in the data is rather 
large. At redshift 2 ^ 3.5 the predicted metallicities are 
much higher than observed, indicating that the combina¬ 
tion of the Equilibrium condition and SHARC assumption 
is invalid at high redshift. 


7.2 Implications of SFR Determined by Halo 
Mass Accretion 

In this section, we look ahead to some important implica¬ 
tions that the toy model has for understanding other cur¬ 
rent modeling techniques. For example, abundance match¬ 
ing, based on galaxy stellar mass and halo mass (or a related 
quantity such as peak circular velocity) leads to the simplest 
models relating galaxies to their host dark matter halos. 
As mentioned in §1, such models predict galaxy correlation 
functions in good agreement with observations both nearby 
and out to high redshifts. But the clustering and galaxy 
content of dark matter halos are known to be a function of 
more than just the mass or circular velocity of the halos|3 
In particular, the formation time and concentration of ha¬ 
los appear to play a major role. Halos of much lower mass 
than the typical mass A4*{z) collapsing at a given redshift 
2 are much more strongly clustered if they formed at high 
redshift than similar-mass hal os that formed at low redshift 
llGao, Springel White|2005l), a phenomenon known as “as- 
semblv bias le.g.. lMo. van den Bosch fc Whit j|2010l'l . Halo 
concentration, Cvir = Rvir/Rs (where Rvir is the virial radius 

^ From the earliest papers on cold dark matter llBlumenthal et 
I1984 lFabeJll984 IPrimacklEgii) it was clear that dark matter 
halos would be characterized by a second parameter beyond mass 
such as overdensity, which is related to formation time. 
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and Rs is the NFW dNavarro. Frenk fc Whitdll997l ~) scale ra¬ 
dius), is related to halo formation time, with halos of higher 
concentration at fixed ma ss forming earlier (|Bullock et alJ 
I2OOII : IWechsler et alj|200^ . and with higher-Cvir halos with 
mass Myir <g M* bei n g much more c l ustere d than average 
llWechsler et al.lboOGl i. iBullock et alJ ll200lll had suggested 
that a natural association is that high-Cvir halos host old, 
red galaxies, and lower Cvir ha los host young, blue gala xies. 
This idea was rediscovered bv iHearin fc WatsonI ll2013l . and 
subsequent papers), who showed that filling halos according 
to this prescription correctly predicts the observed cluster¬ 
ing of red and blue galaxies in the SDSS. It seems surprising 

that such a simple presc ription should work so well. _ 

_ More recently, iHearin, Behroozi fc van den Bos'dil 

1 I 2 OI 5 I) showed that the mass accretion rate of dark matter 
halos at low redshift shows a signal ve ry much like the 
observed two-halo “galaxy conformity” (iKauffmann et alJ 
l2013l : iHearin. Watson fc van den Boschll20l4 ). namely that 
quenched central galaxies tend to lie in quenched regions 
as large as 4 Mpc. A natural explanation for this finding 
would be that the star formation rate in central galaxies is 
closely connected with the mass accretion rate of their host 
halos. 

Figures [7] and [5] show that the sSFR predicted by the 
SHARC assumption, in which the SFR is proportional to the 
host halo mass accretion rate, is in rather good agreement 
with the observations, especially at 2 ~ 0. In more elabo¬ 
rate models of galaxy formation and evolution such as semi- 
analytic models (SAMs) and hydrodynamic simulations, the 
star formation rate, morphology, and other galaxy properties 
are assumed to result from a complex interplay between gas 
inflows and outflows regulated by stellar and AGN feedback 
and other processes, involving much recycled rather than re¬ 
cently accreted gas. If star-forming galaxies at low redshifts 
are indeed nearly in equilibrium, then the SFR will in fact 
be driven by halo mass accretion, which may represent the 
net result of complex processes considered in more detailed 
galaxy formation models. A close connection between halo 
accretion and st ar formation may help to explain the success 
of age matching llHearin fc Watsonll2013f) and the agreement 
between halo MAR conformity and galaxy conform ity obser¬ 
vations (IHearin. Behroozi fc van den Boschll2015ll . 


7.3 Outlook 

Several modifications can add realism to the simplified 
model considered here: 

• Instead of assuming that the stellar-to-halo mass re¬ 
lation (SHMR) is redshift independent, use the evolving 
SHMR implied by abundance matching to connect halo 
MAH to galaxy SFR, using equation ([5]). 

• Instead of assuming that a galaxy is always in equilib¬ 
rium, assume alternatively that the gas mass grows from 
high redshifts down to 2 ~ 4 - i.e., the bathtub fills in the 
early universe. What early universe combinations of galac¬ 
tic gas mass growth and evolving SHMR predict SFRs and 
metallicities in agreement with the rapidly improving obser¬ 
vations? 

• Explore how changing the assumptions regarding gas 
penetration efficiency £eS leads to different dependance on 


halo mass and redshift of mass-loading factors and metallic- 
ity growth. 

• Instead of assuming for simplicity that gas in outgoing 
winds has the ISM metallicity and that freshly accreted gas 
has zero metallicity, compare predictions from modified as¬ 
sumptions with improving data on galactic gas metallicity 
at various redshifts. 

• With the equilibrium condition Mgas.iSM = 0, the gas 
depletion time scale tdepi = Mgas/SFR is just proportional 
to SFR“^, which implies that the slope of the tdeoi to 


pers (ISareent et al. 

I2014 

iGenzel et al. 2015 

) find 


that the slope of the tdepi to sSFR relation is roughly —0.5 
for main sequence galaxies at 2 = 0 to 3. This suggests re¬ 
laxing the equilibrium condition to treat excursions about 
the main sequence. 

• In this paper we considered star formation of central 
galaxies in halos of mass ~ 10^^ to ~ Are there 

simple assumptions that will allow extension of the model 
considered here to more massive galaxies including quench¬ 
ing, and to less massive galaxies including satellites? 


Even without these modifications, we are finding it use¬ 
ful to compare outputs from the simpl e model described 
here to those from a semi-analytic model (IPorter et al.ll2oi3 : 
iBrennan et al.l I 2 OI 5 I I run on the same Bolshoi-Planck ha¬ 
los. We are also comparing the model with zoom-in hydro - 
dynamic galaxy simulations (such as IZolotov et al.l l2015l l . 
We expect that such comparisons will help to improve both 
SAMs and simulations as well as this simple model. 
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